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Asymptotic Behavior of a Class of Stochastic 
Differential Equations 

Zaitang Huang 


Abstract- In this article, we address some conditions on 
invariant measure of Markov semigroups which ensures 
stochastic bifurcations of a wide class of stochastic differential 
equations with fractional Brownian motion. This leads to 
sufficient conditions on Hurst parameter, drift and diffusion 
coefficients for a stochastic bifurcations of the families of 
random dynamical systems.According to the Holder coefficient 
of the diffusion function around the singular point and Hurst 
parameter, we identify different regimes. Interestingly, for the 
first time it is found that the Hurst parameter affects both 
bifurcation conclusions and large deviations which is 
significantly different from the classical Brownian motion 
process. This fact is due to the long-range dependence (LRD) 
property of the fBm. 

Index Terms: Random dynamical systems; Stochastic 

stability; Stochastic bifurcation; Invariant measure; Fractional 
Brownian motion 

I. Introduction 

The classical central limit theorem (CLT) reveals that the 
probability distribution of the sum (or average) of many 
independent and identically distributed (i.i.d.) random 
variables with finite variance approaches the normal 
distribution[l-3].For this reason Gaussian models have been 
widely employed in many fields, and properties of Gaussian 
processes are characterized by second-order statistics, such 
as variance and correlation. In practical applications, 
however, since in many problems related to network traffic 
analysis, mathematical finance, and many other fields the 
processes under study seem empirically to exhibit the 
selfsimilar properties, and the long-range dependent 
properties, and since the fractional Brownian motions are the 
simplest processes of this kind [4- 10]. It is suggested that [11] 
the fractional processes provide a better description of these 
random variables with properties of being infinite variance, 
non-Gaussian or non-stationary. Moreover, it is worth noting 
that the noise term in the stochastic system is used to describe 
the interaction between the (small) system and its (large) 
environment. The non-independence over disjoint time 
intervals in noise term applied by the environment to the 
system makes the fBm more useful since it can exhibit 
long-range dependence. For this reason, and also because the 
fBm includes the fundamental classical Brownian motion as 
a special case when H = 1/2, there has been an increasing 
interest in the research activity related to the fBm itself and 
the stochastic systems driven by it. 
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A large class of phenomena can be described using 
stochastic differential equations(SDE) drive white noise, in 
such distinct domains as economics, physics or biology. 
Many analytical methods have been developed to 
characterize the solutions of such systems. In particular, the 
widely applied theory of stochastic bifurcations [12] allows a 
qualitative characterization of the asymptotic regimes of 
random dynamical systems. This theory is very handy to 
study nonlinear stochastic differential equations and is used 
to characterize the asymptotic behavior of complicated 
systems. Investigating the impact of fractional Brownian 
motion noisy perturbations on such systems is hence of great 
interest, and is currently an active field of research. In recent 
years, it has been shown in many different areas that applying 
fractional noise on a random system can lead to many 
counter-intuitive phenomena, such as fractional 
noise-induced stabilization [13] and control[14]. From a 
mathematical perspective, understanding the interplays 
between white noise and fractional noisy perturbations is a 
great challenge, with many applications. Several tools have 
been introduced, ranging from the theory of random 
dynamical systems (RDS) [15,16], the study of moment 
equations [17,18] or multiscale stochastic methods for 
slow-fast systems [19]. Unfortunately no systematic method 
in the flavor of bifurcation theory for the analysis of the 
dynamics of nonlinear SDE drive fractional white noise exist, 
and this is the central problematic of the present manuscript. 

We focus here on the dynamics of random dynamical 
systems induced parametrized one dimensional stochastic 
differential equation driven by the fractional Brownian 
motion(fBm). The question we address is how the interplay 
between the Hurst parameter, drift and the shape of the 
diffusion function affects the behavior of the system. An 
important contribution of RDS theory to the field of 
stochastic bifurcations is to distinguish between 
phenomenological(P) bifurcations and dynamical (D) 
bifurcations. One is the phenomenological approach favored 
by physicists and engineers based on the qualitative changes 
of stationary measured. e., the stationary probability density 
of the response. The other is the dynamical approach favored 
by mathematicians based on the qualitative changes of 
stability of invariant measures and occurrence of new 
invariant measures for random dynamical systems. Each 
approach has its advantages and the two approaches can be 
regarded as complementary to each other [12,20]. We 
investigate the questions of stochastic stability and 
bifurcations)? and D) in a comprehensive manner, using 
combinations of appropriate tools such as invariant measure, 
Lyapunov exponent, stationary measure and Lyapunov 
functionals. Mathematically, the choice of a vanishing 
diffusion coefficient provides a framework to study subtle 
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competitions between the drift and the fractional noise, 
leading to a rich and generic phenomenology. 

The paper is organized of as follows: In Section 2, we give 
some Definitions. We study the stochastic bifurcation in 
Section 3 and the two examples are given to illustrative our 
presented results in Section 4. 

III. PRELIMINARIES 

In the section, we present some preliminary results to be 
used in a subsequent section to establish the stochastic 
stability and stochastic bifurcation. Before proving the main 
theorem we give some lemmas and definitions. 

In this paper, we analyze the behavior of the solutions for a 
class of stochastic equation with fractional Brownian 
motion(fBm): 

dx l =f(x t )dt + g(x t )dB f, (2.1) 

where is a fractional Brownian motion defined on a 
complete probability space (Cl, F ,P) endowed with the 

natural filtration (F ,)\ °f B ^ .Notice that several classical 

normal forms can be reformulated in this setting and our 
general results will be applied to the cases of the Transcritical, 
pitchfork or Hopf bifurcations in section 3. We assume that 

there exists a singular poin x such that fix ) = g(x ) = 0 . 

The system has a trivial solution distributed as a Dirac 

* 

measure localized at the singular point X , and we shall 
address the stochastic bifurcation of these solutions. In 
contrast to the case of ODEs there are at least five different 
notions of stochastic bifurcation [12] 

Definition 1 (D-Bifurcation) Dynamical bifurcation 
is concerned with a family random dynamical systems which 
is differential and has invariant measure jU a . If there exist a 

constant Oi D satisfying in any neighborhood of Oi D , there 
exist another constant OC and the corresponding invariant 
measure V a ^ p a satisfying V a — > p a as Oi — > Oi D . 

Then the constant Oi D is a point of dynamical bifurcation. 

Definition 2 (P-Bifurcation) Phenomenological 

bifurcation is concerned with the change in the shape of 
density (stationary probability density) of a family random 
dynamical systems as the change of the parameter. If there 
exist a constant CC 0 satisfying in any neighborhood of 0i D , 

there exist other two constant d x , CC 2 and their corresponding 
stationary density function p a ^ , p a satisfying p a and p a 

are not equivalent. Then the constant Oi 0 is a point of 

phenomenological bifurcation. 

Definition 3 (Stochastic Transcritical Bifurcation) 

A family of RDS (< 'P a ) a eS ; on undergoes a Stochastic 

Transcritical Bifurcation at Oi = 0 , if 

(i) for Oi < 0 , (p a has exactly two ergodic invariant 

measure in I {(p a ) : 5 {] which is stable, and fl a = 8 a with a 
random variable a a < OP as - which is unstable, and 


a a — > 0 in probability as (X T 0 , 

(ii) for Oi — 0 , S 0 is only invariant measure and the 
Lyapunov exponent of (p Q with respect to 8 0 vanishes, 

(iii) for a >0 ,(p a has exactly two ergodic invariant 
measure in / ((p a ) : S f) which is unstable, and jl a =8 U 
with a random variable Cl a > OP , which is stable, and 

Cl a — > 0 in probability as (/ -i 0 , 

Definition 4 (Stochastic Pitchfork Bifurcation) A 

family of RDS (<p ) aG2lt on undergoes a Stochastic 
Pitchfork Bifurcation at a = 0 , if 

(i) for d < 0 , 8 0 is only invariant measure of (p a , which 
is stable tor (X < 0 and the Lyapunov exponent of (p {) with 
respect to 8 Q vanishes; 

(ii) for Oi > 0 the system possesses besides 8 0 , which is 
unstable, exactly two more ergodic invariant measure 
f/ a , lC a in I {(p a ) , described by fj! a - 8 , , i = 1, 2, with 

a a 

random variable Cl' a > 0 , cC u < OP a v and p' a , i = 1, 2, 
are stable; 

(iii) We have a' a — > 0 in probability as a X 0, i = 1, 2. 

Definition 5 (Stochastic Hopf Bifurcation) 

In the viewpoint of phenomenological bifurcation: 

The stationary solution of the FPK equation which is 
corresponded to the stochastic differential equation changes 
from one peak into crater. 

In the viewpoint of dynamical bifurcation: 

(i) If one of the invariant measures of the stochastic 
differential equation loses its stability and becomes unstable 
(i.e. two Lyapunov exponent are positive), moreover the 
rotation number is not zero. Meanwhile there at least appear 
one new stable invariant measure. 

(ii) The global attractors of the stochastic differential 
equation change from a single point set into a random 
topological disk (the closure of the unstable manifold of the 
unstable invariant measure). 

If the stochastic bifurcation of a stochastic differential 
equation has the above characters, then the stochastic 
differential equation admits stochastic Hopf bifurcation. 

Definition 6 (Stochastically Stable) The trivial solution 

x(t,t Q ,X 0 ) of stochastic differential equation is said to be 
stochastically stable or stable in probability if for every pair 
of S £ (0, 1) and a > 0 , there exists a 5 = 5{e , 8, t 0 ) > 0 
such that 

P{\x(t , t 0 , X () )| < a for all t > t 0 } > 1 - £ , 

whenever |x 0 1 < 8 . Otherwise, it is said to be stochastically 
unstable. 

III. STOCHASTIC STBILITY AND BIFURCATION 

Now that we characterized in detail the solutions of 
equation (2.1), we are in a position address the dynamics of 
three prominent bifurcations occurring in the study of 
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random dynamical systems: the transcritical, pitchfork and 
Hopf bifurcations with power diffusion coefficients. 

A. Stochastic Transcritical Bifurcation 

In this section we first study the dynamics of an SDE with 
fractional Brownian motion(fBm) given by the normal form 
of the transcritical bifurcation, before addressing the 
universality of these behaviors. 

We consider the stochastic differential equation with fBm 

dx t = ( ax t -x 2 )dt + cjfdB" 

( 3 . 1 ) 

Theorem 3.1 For k = 1, let us denote by x () the 
stationary solution. We have: 

( 1 ) When 0 < H < 1/2 and: 

(1) For CC < 0 , the fixed point X {) of equation (3.1) is stable 

in probability and the two invariant measures are both 
measurable; 

(ii) For CC > 0 , the fixed point x 0 of equation (3.1) is 
unstable in probability and the invariant measure is both 
F measurable; 

( 2 ) When H = 1 / 2 , and: 

cr 2 

(i) For CC — — < 0 , the fi xed point X Q of equation (3.1) 

is stable in probability and the two invariant measures are 

both F 0 measurable; 

—00 ’ 

a i 

(ii) For CC — — > 0 , the fi xed point X 0 of equation (3.1) 

is unstable in probability and the invariant measure is both 
measurable; 

(3) When 1/2 < H < 1 , and (X E , the fi xed point 

X 0 of equation (3.1) is almost surely exponentially stable. 

Moreover, any solution almost surely reaches zero in finite 
time. 


Proof. We now calculate the Lyapunov exponent for 
each of these measures of the SDE(3.1). The linearized RDS 
V t = D<p(t,C0,x)v satisfies the linearized SDE(3.1) with 
fBm 


dv t = [cc-<j Ht 2H 1 - 2(cp{t, co)x)\v l dt + a\\dB ^ , 

hence 


D(p(t,oo,x)v = vexp 



+ oBf (co) 


A( . =lim-log||D (t,co,x)v\ 

?-» 00 f M Y II 

2 i 

= a~ — t 2H ^ -21im- [ (< p(t,co)x) N l ds 

^ t — >oo ^ JO 

2 

= a limr H1 -2Ex,x 5 

2 »-*» 

provided the IC x 0 eL'(P) is satisfied. 


From SDE(3.1) with fBm, we can obtain the solution of 

SDE(3. 1) which is explicitly given by which is solved by 

(p a {t,C0)X = 


xe 


at-dl t 2H +cjS « (*,) 
2 ' 


1 + xf 

J( 


t as-—s 2H +oB f (®) 

e 2 

0 


ds 


We now determine the domain D a (t,C0) and the range 

R a (t,C0) of (p a {t , 0)) : D a (t,oj) R a (t,co) 


We have 


D(t,cd) = 


(-d a ( t , co), 00 ), t > 0 , 

7 = 0 , 

(- 00 ,d(t,co)), t< 0 , 


where 


dJt,(o) = 


f 


I, 


t as-—s 2H +crB 2l {co) 

e 2 

0 


ds 


>0, 


) 


and 

R a ( t , co) = D a (-t, O(t)co) 


(-°°, r a (0 «))» 

(~r(t,co), 00 ), 


t > 0 , 
7 = 0 , 

7 > 0 , 


—2 Jo ((pis, oj)x)dsj. where 

Thus, if jJ„, = S Xo(a) is a ^-invariant measure its Lyapunov 
exponent is 
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r (t, CO) = d r/ i-t, Git)co ) 


at—F—t 2H +CTB" 


f 

2 

A 

rt (as—F^ s 2H +crBd{co)) 



e 2 ds 



Jo 


V 


) 


> 0 . 


We can now determine 
(i) For 0<H < 1/2, 


[ 0 ,d a (co)], a> 0 , 
0 , a = 0 , 
i-dUco), 0 ], «< 0 , 


(ii) For H =1/2, 


E =< 

CC 


[O ,d„ia>y\, oc 


cr 


O, 


CT 


CX = 


(— cC(cfc>),0], « < 


CT 


(iii) For 1/2 < H < 1, 

E a = (“C (®)> 0], a, cr e 91, 

where 

0<d*(o) = 


1 


CT" 


f 

CT 2 o 

A 


l*±oo (as s 2H +aB^ (co)) 



\ e 2 ds 

JO 


V 


J 


< oo. 


(l)For0 < H < 1/2 and: 

(i) When G; e 91 , the inclusion closed(IC) /// w = S 0 is 
trivially satisfied and we obtain 

A(ju“ ) = CC 

so Hi i s stable for CC < 0 and unstable for CC > 0 . 


Fokker-Planck-Kolmogorov(FPK) equation[4] 


d P a (t x ) d d,ax-—x-x)p a (t,x)) 


dt dx 

d 2 ((aJ x(f)(s, t)ds ) p a (t, a)) 


5a 2 


-0 


which has the unique probability density solution 
P a (t, A) [4], 

Since 

/•GO 

E a x = E(d") = [ xp a (x)ds<cc, 

th JO 

the IC is satisfi ed. The calculation of the Laypunov 
exponent is accomplished by observing that 


d a _(0 t 6)) = 


2 

{at — —t 2H +(jB? {co)) 


y/\t) 


I 


t {as-E S 2H {a))ds y/{t) 

(e 2 


where 


v(t) = J 

J - 


t {as~~-,$ 2H +<jB^ (co)) 


ds. 


Hence by the ergodic theorem 

E(d a ) = lim-log^(f) = a, 
t 

finally 

A(ju%) = -a < 0. 

(iii) When a < 0 , jLl“ 0J = 8 d „ (o)) is F measurable. 
Since L (d d )—L (d a j and N is even, 

E(-d*) = -E(d: a ) = a 


thus 


A (// 3 " ) = -a > 0. 


(ii) 

When 

55 

V 

O 

Kco = 8 c i" (a) is F -l measurable. 

so /.// 

is 

hence 

the 

density 

p a of p a = E „ satisfies the 

ft 

a 2 
a 

>0. 


Hence, we have a D-bifurcation of the trivial reference 

2 

cr a 

measure o Q at OC D = and a P-bifurcation of the V ± at 

dtp. Then the SDE(3.1) with fBm admits stochastic 
Transcritical bifurcation. 

(2) H = 1/2 and: 

(i) When CC, CT £ 91, the inclusion closed(IC) for 
JLl\ f0 = <b' (l is trivially satisfied and we obtain 

a(h“) = a ~~Y 

^2 

ble for a < 0 and unstable for 

2 
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(ii)When a - — > 0 , a 

2 


5 d“(0» im- 


measurable. 


hence the density p" of p a = E//“ satisfies the 
Fokker-Planck-Kolmogorov(FPK) equation[4] 
a 2 

dp° d«ax- — x-x-)p°) g 2 (x 2 p « } 

dt dx dx 2 

which has the unique probability density solution 

W -i e xp(-% x>0, 


p = N„x 

where 


(3) 1/2 < H < 1, the inclusion closed(IC) for /u" m = S Q 
is trivially satisfied and we obtain 

) = - 00 , 

so /// , fl“ and ju“ are stable for a e S .H. 

Let us now further describe the dynamics of the solutions 
= 0 (3.2) as a function of H. 

U Case H — 1/2 : We observe that for any CC> 0, the 
unstable deterministic fixed point becomes asymptotically 
exponentially stable when the noise parameter (7 is large 


<J 


ac 1 =r 


(2a) 

V 1 

2 

\o- J 

U J 


a 


enough. For a > , two symmetrical stationary 

2 

cr 2 

distributions appear. For — < a < ( 7“ the stationary 

2 


Since 

poo 

x = E (d “ ) = xp a (x)dx < oo ; 


the IC is satisfied. The calculation of the Laypunov exponent 
is accomplished by observing that 


d a _{.0,co) = 


(at-2—t 2H +aBf 1 (w)) 


v\t) 


rt (as-—s 2H +aB^(a>))ds y/(t) 

(e 2 

J — oo 


where 


HO = j 


f (as-—s 2I, +aB^(a>)) 

e 2 


ds. 


Hence by the ergodic theorem 


finally 


E (d a ) = lim - log y/(t ) = a - — 
r^oo t 2 




2 

(iii) When a-°Y< 0, = d_ d , /(oj . is F_° 

measurable. 

Since L = L a ^, 

2 

E (-d * ) = -E(d: a ) = a - — 

thus 


A(^“) = -(a-y)> 0 . 

Hence, we have a D-bifurcation of the trivial reference 
measure § 0 at a D = — and a P-bifurcation of the V at 

c p 

ap — . Then the SDE(3.1) with fBm admits stochastic 

2 

Transcritical bifurcation. 


distribution concentrates at zero and has a non increasing 
density diverging at zero. For a><7~, the probability 
density function of the stationary distribution vanishes at zero 
and has a unique maximum reached for x ± \la — fJ 2 . There 
is hence a qualitative transition at a = cr 2 , or P-bifurcation. 
In comparison with the deterministic bifurcation, the loss of 
stability is delayed and noise tends to stabilize the saddle 
point. 

• Case 0 < H <1/2: the deterministic picture is 
qualitatively and quantitatively recovered: for a< 0, X 0 is 
stable in probability and is the unique stationary solution, and 
for a > 0 is unstable in probability, two additional 
stationary solutions appear. 

• Case 1/2 < H < 1 : the solution X Q is always stable in 
probability. This result can appear relatively surprising at first 
sight. Indeed, in the deterministic case, (/ is the exponential 
rate of divergence from the solution 0. However, adding a 
(possibly small) diffusion term proportional to |x| dB H (t) 

with 1/2 < H < 1 stabilizes X 0 in probability whatever the 

value of the noise intensity. This observation, added to the 
fact that the solution is not exponentially asymptotically 
stable (though stable in probability) raises the question of 
how this convergence occurs. Starting from a positive initial 
condition, we observe that the solution is evolving in the 
half-plane x > 0 and does not show any extinction clue. 
However, it suddenly reaches zero where it is absorbed after 
that random transient phase. This perfectly illustrates a 
typical behavior of the solutions of the transcritical equation 
for 1/2 < H < 1 : the extinction time was shown to be almost 
surely finite. Moreover, trajectories that did not hit zero are 
distributed according to a quasi -stationary distributions as 
long as 1/2 < H < 1 . Interestingly, for the first time it is 
found that the Hurst parameter affects both bifurcation 
conclusions and large deviations which is significantly 
different from the classical Brownian motion process. This 
fact is due to the long-range dependence (LRD) property of 
the fBm. 
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Theorem 3.2 For k < 1, 0 < H < 1 , then the fixed point 
X 0 of equation (3.1) is asymptotically stable in probability 

whatever and (7 X 0 . 

Proof. Case 0 < H < 0.5 : Let us consider V(x) = x k 
as a Lyapunov functional for the dynamics. This function is 
clearly C\ (91) , i.e. it is twice continuously differentiable 
except at X — 0 .Moreover, we have: 

LV < kx kl (ax—x 2 ) — k(l — k)a 2 Ht 2H1 x 3k2 . 
Since we assume that 0<//<0.5and —k(l — k)Ht~ n 1 

3k-2 k + 1 

X at 0, the leading term close to 0 is of order —X 
which is strictly negative for sufficiently small x. More 

precisely, there exists r such that V EC^(U r ) where U r 
is the open ball of radius r and such that LV < 0 for all 
X £ U r . Moreover, for any £ such that 0 < £ < r we have 
V (r) and LV < —C s < 0 for all r>X>£. Theorem of 
[21] therefore applies and concludes the proof of the stability 
in probability of the solution X 0 whatever the parameters 
«£ s . R, <x > 0 ■ 

Case 0.5 <H <1: Let us consider V(x) = x as a 
Lyapunov functional for the dynamics. This function is 
clearly C, (9?) , i.e. it is twice continuously differentiable 
except at X = 0 . Moreover, we have: 

LV < kx k (ax-x 2 )-k{ 1 - k)a 2 Ht 2Hl x 3k2 . 

Since we assume that 0 < H < 0.5 and —kx k+l at 0, the 
leading term close to 0 is of order 

—k(l—k)(J~Ht~ H 1 X 31 " which is strictly negative for 
sufficiently small X . More precisely, there exists r such 
that V £ C, (U ) where U r is the open ball of radius 
r and such thatL V < 0 for all X £ U r . Moreover, for any 
£ such that 0 < £ < r we have V (r) and 
LV <—C E < 0 for all r>X>£ . Theorem of [21] 
therefore applies and concludes the proof of the stability in 
probability of the solution X 0 whatever the parameters 
CC £ 9?, <7 > 0 . When H =0.5, then k £ [1 / 2, 1) . 
Theorem 3.3 For k > 1, let us denote by X 0 the 

stationary solution X =0 a.s. for all t . We have: 

(l)For 0 < H < 1/2 and: 

(1) When (X < 0 . the fixed point X 0 of equation (3.1) is 
stable in probability and the two invariant measures are both 
F measurable; 

(ii) When CC >0 , the fixed point X 0 of equation (3.1) is 
unstable in probability and the invariant measure is both 
F {) ' measurable; 

(2) For 1/2 < H <1 and Ol £ 91, the fixed point X 0 of 
equation (3.1) is not stable in probability. 


Proof. Case 0 < H < 1/2 : Let us consider 

V (x) = x 2 as a Lyapunov functional for the dynamics. This 
function is clearly C\ (31) , i.e. it is twice continuously 
differentiable except at X = 0. Moreover, we have: 

LV < 2x(ccx - x 2 ) + a 2 Ht 2H - l x 2k . 

Since we assume that 0 < H <1/2 and — 2x at 0, the 

term LV is equivalent to 2 OCX" and hence locally has the 
sign of OC. Force < 0, we can directly apply theorem [21]. 
For CC >0 , we use V (x) = — log(x) again and conclude 
that 0 is not stable in probability. 

Case 0.5 < FI < 1 : Let us consider V (x) = X 2 as a 
Lyapunov functional for the dynamics. This function is 
clearly C/ (31), i.e. it is twice continuously differentiable 
except at A =0. Moreover, we have: 

LV < 2 x(ax - x 2 ) + <j 2 Ht 2HX x k . 

Since we assume that 1 / 2 < H <1 and — 2x at 0, the term 
LV is equivalent to c 2 Ht 1Hl x 4 and hence LV > 0 , we 
can directly apply theorem [21], For a £ 31 we conclude that 
0 is not stable in probability. 

Ii. Stochastic Pitchfork Bifurcation 

In this section we first study the dynamics of an SDE with 
fractional Brownian motion(fBm) given by the normal form 
of the pitchfork bifurcation, before addressing the 
universality of these behaviors. 

We consider the supercritical stochastic pitchfork 
bifurcation with fBm 

dx, = (ax, - x, 3 )dt + a | x, |* dB kI (3.3) 

Theorem 3.4 For H = 1/2 , let us denote by X 0 the 
stationary solution X = 0 a.s. for all t. We have: 

(l)For k = 1 and: 

0.2 

(1) When a < 0 , the fixed point X 0 of equation (3.1) 

is stable in probability and the only invariant measure is 

CT 2 

(ii)When a > 0 , the fixed point X 0 of equation (3.1) 

is stable in probability and the three forward Markov measure 

Ml = and v tco = 5 d a («) . where 

1 

dfm):=( 2j° 

(2) For k< 1, the fixed point X 0 of equation (3.1) is 
asymptotically stable in probability whatever CC £ 91 and 
a ^0 . 

(3) For A: > 1, x 0 is 

-stable in probability for CC <0 
-unstable in probability for CC > 0 

Proof. (l)Casek=l. 
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We now calculate the Lyapunov exponent for each of 
these measures of the SDE(3.3). The linearized RDS _ jy ^ ^ O(t)O)) 

V t = D (t, 0), X)V satisfies the linearized SDE(3.3) with 
ffim 

dv t =(a- cr 2 Ht 2H 1 - 3 (cp(t, a>)x) 2 v t dt + crv t dB" , (3.4) 

hence 

f _2 

D(p(t,o), x)v- vexp 


(-°°» r a o))» 

= \ m, 

(- r(t , co), go). 


t > 0, 
7 = 0 , 
t > 0, 


at-—t 2H + (tB' 1 (co) 


where 


V 


-3 ^((p(s,co)x)) 1 ds 


r a {t,co) = d a [-t,e(t)co) 


Thus, if JU (0 = S' is a cp -invariant measure its Lyapunov 
exponent is 

X(ju) = lim - log \\D<p(t, co, x)v|| 


e 


O fj 

at t+oB, 

2 ' 


f 


2 

2 

<J 


- 2 lim co)x) N 1 )ds 


1 , 


t . as s+crBf (cd) 

e 


ds 


xl/2 


> 0 . 


= a 3 Ex, 

2 0 

provided the IC x 0 e L 1 (P) is satisfied. 

From SDE(3.1) with fBm, we can obtain the solution of 
SDE(3.1) which is explicitly given by which is solved by 


We can now determine 

0,d- (co) 


E„ ={ 


(p(t,(D)X = 


xe 


at-—t+<jB?(co) 
2 ' 


l + 2rf 

Ji 


2 y/2 

O' r,H , 


0, 


a 


a > 

a 2 

T 


cr 


(~d + a (co), 0], a< ~Y’ 


2 <•/ 2(as s+aB s (co)) 


ds 


V 


where 

0<d;(w) 


We now determine the domain D a (t,CO) and the range 

R a (t,co) of cp a (t , CO) : D a {t,CO)^> R a (t,CO ) 

We have 

D a it, co) 


where 


(~d a it,co),oc). 

t > O 

m, t = o. 


( go, d tx it, co)). 

t < o 

1 



f 

(a 2 1 





2 

f" 2 } ds 

Jo 


V 


) 


y/2 


< 00 . 


From Proposition 2.1, there are ergodic invariant measures 
of SDE(3.1) with fBm as follows: 

2 

(i) For a < 0, the only invariant measure 

2 




cr 


f 

cr 2 H 

\ 


rt (as s+crB) 1 ( a))ds 


2 

\e 2 

Jo 


V 


) 


(ii) For CC — — > 0, the fixed point X 0 of equation (3.1) 

is stable in probability and the three forward Markov 
measure //“ = S 0 and v“ a = S d ^ a) , where 

1 


> 0 , 


and 


djoc) := 2j 


0 


2at-a^t+aB^ (co) 


\ 


dt 
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We have E d~(co^ = a — — . Solving the forward 


Fokker-Planck-Kolmogorov(FPK) 

a 

L*p a = - 


equation[12] 


f 

r 

2 

k 

\ 



a 

3 



ax- 

— x- 

X 

p 



2 



vv 


) 

J 


8x 

yields 

(i) p a = S 0 for all a, a, 

2 

a 

(ii) for a > 0 

2 


+ 


d\(x z p a ) 

dx 2 


= 0 


da 


2 a 

N„x^~ 


c 


exp 


2 a 


\ 


"-ZT* 
V cr 


x > 0, 


C ia{ X ) = ( la{- X )- 

Naturally the invariant measure = 8 d are those 

corresponding to the stationary measures q±. Hence all 
invariant measure are Markov measures. 


The two families of densities 


(«*), 


a-— > 0 
2 


clearly undergo 


a P-bifurcation at the parameter value ap — — 

2 

By Proposition 2.1, the Lyapunov exponent of the 
linearized SDE(3.3) with respect to the three measure is 

_2 

(i) for JLl a = S 0 ,/ t ( fJL a j = CC 


a 

~2 


(ii)for K 0 > = 8 , 


d a (ty ) 9 


A(v a ) = a- — -3Ed 2 a =-2cc + <j 2 . 

Hence, we have a D-bifurcation of the trivial reference 
a 2 

measure 8 0 at a D = and a P-bifurcation of the v" at 

g .2 

a,, = .Then the SDE(3.4) with fBm admits stochastic 

2 

pitchfork bifurcation. 

(2)Case k< 1. 

Let us consider V (x) = x k as a Lyapunov functional for the 

dynamics. This function is clearly C° (9?), i.e. it is twice 

continuously differentiable except at X = 0 . Moreover, we 
have: 


LV <kx k -\ax-x i )~ 


k(l-k) 
2 


-a x 


2 2k— 2 


Since we assume that k < 1 and —kx + at 0,the leading term 
close to 0 is of order —k(\ — k)<J 2 X 3k 2 / 2 which is strictly 
negative for sufficiently small X. More precisely, there 
exists r > 0 such that V e C" (iJ r ) where U r is the open ball 
of radius r and such that LV < 0 for all X £ U r .Moreover, 
for any 8 such that ()<£■< r we have V (r) and LV < —C E 


N. 
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< 0 for all r > X > £ .Theorem of [21] therefore applies and 
concludes the proof of the stablility in probability of the 
solution X 0 whatever the parameters CC £ ‘J? , (7 > 0 

and A' G [1/2,1). 

(3)Case k> 1 

Let us consider V(x) — x 2 as a Lyapunov functional for the 

dynamics. This function is clearly C/ (91 ) , i.e. it is twice 
continuously differentiable except at X = 0 . Moreover, we 
have: 

2 

LV < 2x[ax — x 2 ^ + ^- x 2k . 

Since we assume that k > 1 and — 2x 4 at 0,the term L V is 
equivalent to 2 QOT and hence locally has the sign 
of a. For CC < 0, we can directly apply theorem of [21], 
For CC > 0 ,we use V(x) — — log(x) again and conclude that 0 
is not stable in probability. 


Theorem 3.5 For 0 <H ^ 1 j 2 , let us denote by Xq the 
stationary solution X = 0 a.s. for all t. We have: 

(l)For k = 1 and: 

(1) CC < 0, the fixed point X 0 of equation (3.1) is stable in 
probability and the only invariant measure is p ' n = S Q ; 

(ii) Of > 0, the fixed point X 0 of equation (3.1) is stable in 
probability and the three forward Markov measure 
and//“ = S () and v“ w = S da(a) , where 

(2) For k < 1 and: 

(i) CC < 0, X Q is stable in probability and no other 
stationary solution exist. 

(ii) CC > 0, Xq is unstable in probability for a > 0 and 
stationary solution exist. 

(3) For k > 1 and: 

(i) CC < 0, X 0 is stable in probability and no other 
stationary solution exist. 

(ii) CC > 0, Xq is unstable in probability for CC >0 and 
stationary solution exist. 


Proof. (1) Case k=l 

From SDE(3.3) with fBm, we have the linearized RDS V t 
= D(p(t, CO, x) V satisfies the linearized SDE(3.3) with fBm 

dv t ={^a-<j 2 Ht 2H X -3[cj){t,(o)xy^v t dt 

+ av t dBf , 

hence 

D(j)(t,C 0 ,x)v = vexp 

( 2 f ' 

at t 2H +<jBf (<y)-3| ((z>(v,&>)x)~ els 

v ^ 0 


Thus, if p co = 8 x is a (p -invariant measure its Lyapunov 

exponent is 
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A(//) = lim-log||Z)^(f,ty,x)v|| 

= a-^-t 2H ~ l -21im-| o (^(f,®)x) iV 1 ds 

= a-—t 2HX - 3Ex 0 2 
2 0 

provided the IC X 0 E ll ( P ) is satisfied. 

From SDE (3.6) with fBm, we can obtain the solution of 
SDE(3.6) which is explicitly given by which is solved by 

(p a (t,co)x = 


xe 


at~t 2H - l +*B?(a>) 



v y 

From Proposition 2.1, there are ergodic invariant measures 
of SDE (3.6) with fBm as follows: 

(i) For a < 0, the only invariant measure is //“ = S () ; 


(ii) For a > 0, the fixed point X Q of equation (3.1) is 
stable in probability and the three forward Markov 
measure /// = (),, and v“,, = S, , , , where 

‘ CO U ±,CO d a (G)) ’ 


d « M 2 Jl e 


2 " * +crB H ( co) 


dt 


We have Er// (<y) — CC — . Solving the forward 

Fokker-Planck-Kolmogorov(FPK) equation[4] 


(t,x) ^[ax-x 3 )p a (t,xfj 


dt 


dx 


d 2 (x 2 a 2 Ht 2H ~ l p a (t,x)) 

d? 

yields 

(i) p a = <5 0 for all a, a. 


= 0 


(ii) for a > 0 q a (x) = q + a (-*) , x > 0. 

Naturally the invariant measure V+ ca = S d f 0J j are those 

corresponding to the stationary measures q*. Hence all 
invariant measure are Markov measures. 

The two families of densities {q ' u ) // >0 clearly undergo a 
P-bifurcation. 

By Proposition 2.1, the Lyapunov exponent of the 
linearized SDE (3.6) with respect to the three measure is 

(i) for ju a — S 0 ,h{/J a ^)= a , 


N. 
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(ii) for v a ±m = S dJa) (y a ) = a- 3E d 2 a = -2a. 
Hence, we have a D-bifurcation of the trivial reference 
measure S fj at a D = 0 and and a P-bifurcation of the v" at 

2 

a = — Then the SDE (3.4) with fBm admits stochastic 

p 2 

pitchfork bifurcation. 

(2) Case k< 1 

Let us consider V(x) = x k as a Lyapunov functional for 

the dynamics. This function is clearly C° (93) , i.e .it is twice 

continuously differentiable except at X = 0 . Moreover, we 
have: 

L V < kx k - x (ax - x 3 ) - k (l - k ) a 2 Ht 2H - x x u ~ 2 . 

Since we assume that k < 1 and — 1 

X 3k ~ at 0, the leading term close to 0 is of order kox which 
is strictly negative for sufficiently small X and a < 0 ;the 
leading term close to 0 is of order k(XX which is strictly 
positive for sufficiently small X and a > 0, we can directly 

apply theorem of [2]. For a < 0 , X f) is stable in probability; 
for a > 0 , X 0 is stable in probability. 

(3) Case k> 1 

Let us consider V(x) = x 2 as a Lyapunov functional for the 
dynamics. This function is clearly C° (93) , i.e. it is twice 
continuously differentiable except at X = 0 . Moreover, we 
have: 

2x(ax - x 3 ) + 2 a 2 Ht 2H - l x 2k . 

Since we assume that k > 1 and 2 <j 2 Ht 2H ~ l X 2k at 0 , the 
term L V is equivalent to 2ax 2 and hence locally has the 
sign of a . For a < 0 ,we can directly apply theorem of [21], 
X 0 is stable in probability. For a > 0 , we use V(x) = -log(x) 
again and conclude that X 0 is not stable in probability. 


Theorem 3.6 For 1/2 < H < 1, , let us denote by X 0 
stationary solution X = 0 a.s.for all t. We have: 

(1) For k < 1, X = 0 is asymptotically stable in probability 
whatever a E 93 and (j X 0 . 

(2) For k > l,X = 0 is unstable in probability 
for a e 93 and cr ^ 0 . 

Proof. (1) Case k= 1 

From SDE (3.3) with fBm, The linearized 

RDS V t = D ip ([,OJ,x)v 


satisfies the linearized SDE(3.3) with fBm 

dv t =(a-cr 2 Ht 2H l -3 ((p(f,®)x)" Wrfr 


(3.6) 


+ (JvdB n , 


hence 
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D(p(t,oo,x)v = vexp 

( 2 t \ 

at-—t 2H + oBf 1 (co)- 3j {(p{s,co)xf ds 
V 2 0 

Thus, if /l a = 8 k is a <p -invariant measure its Lyapunov 
exponent is 

Z(//) = lim-log|| J D^(f,ty,x)v|| 

= a-^-t 2H ^ -21im-| o (^(f,r9)x) iV 1 ds 

= a-—t 2H ~ 1 - 3Ex 2 

2 0 

Since we assume that k > 1 , then 2,(//) < 0 for CC £ )R . 

(2) Case k< 1 

Let us consider V (x) = x k as a Lyapunov functional for 
the dynamics. This function is clearly C' l (iR) , i.e. it is twice 
continuously dierentiable except at X = 0 . Moreover, we 
have: 

L V < kx k ~ l (ax-x 3 )-k(l-k) a 2 Ht 2H - l x 3k ~ 2 . 

Since we assume that k < 1 and —kx at 0,the leading term 
close to 0 is of order —k(\ — k)<J~X 3k 2 / 2 which is strictly 
negative for sufficiently small X. More precisely, there 
exists r > 0 such that V £ C° (C/ r ) where U r is the open ball 
of radius r and such thatL V < 0 for all X £ U r .Moreover, 
for any £ such that 0 < £ < r we have V (r) and L V < —C E 

< 0 for all r > X > £ .Theorem of [21] therefore applies and 
concludes the proof of the stablility in probability of the 
solution X Q whatever the parameters CC £ ‘J? , <J > 0 and 

k e [1/2,1). 

(3) Case k> 1 

Let us consider V(x) — x 2 as a Lyapunov functional for the 
dynamics. This function is clearly C° (iR) , i.e. it is twice 
continuously differentiable except at X = 0 . Moreover, we 
have: 

2x(ax-x 3 ) + 2a 1 Hr H - 1 x 2k . 

Since we assume that k < 1 and —kx at 0,the leading term 

close to 0 is of order 2 <J~Ht~ H l X~ k which is strictly 
negative for sufficiently small X. More precisely, there 
exists r > 0 such that V e C° (U r ) where U r is the open ball 
of radius r and such that L V > 0 for all X £ U r .Moreover, 
for any £ such that 0 < £ < r we have V (r) and 
LV >c E > 

0 for all r > X > £ .Theorem of [2] therefore applies and 
concludes the proof of the stablility in probability of the 
solution X 0 whatever the parameters CC £ V .R , (7 >0. 

C. Stochastic Hopf Bifurcation 


We study in this section the stochastic Hopf bifurcation 
with multiplicative noise. For simplicity, we will consider 
that the stochastic perturbations are driven by a single with 
fractional Brownian motion. 

We consider the stochastic Hopf normal form with fBm 

dX, = (aX,-Y t +£X t [x 2 +Y 2 )dt} 

+ (cr 1 X t -cJ 2 Y t )\Z t \ k - ld B?, 

1 / (3 ' 7) 
dY t =(aY t + X t +sY t (X 2 +Y 2 )dt J 

+ (cr 2 X t +a l Y t )\Z r \ k - 1 dB? , 

where a and £ correspond to the parameters of the Hopf 
bifurcation, Z t = X t + iY t and a t and ay are two 

parameters globally governing the amplitude of the stochastic 
perturbation. 

The deterministic process Z ( corresponding to X t = 0, 

Y t = 0 for all t > 0 is solution of ( 3 . 7 ) and is univocally 

defined by the fact that the modulus of Z t is null. We denote 

by R t the modulus of Z ( and by R f its argument 

Lemma 1. The modulus of the variable R t = |Z ( | and the 

argument 9 t satisfy the equations: 

dR t = ( aR, + a 2 R 2k ~ l Ht 2H ~ l + cR , 3 ) dt + cr^'dB' 1 

< ( 3 . 8 ) 

d0,={ 1 - 2a l a 1 R 2k ~ 2 Ht 1H ~ l ) dt + o^'dB" 
where R t — \Z t | = yj X~ + Y~ , 0 t given by the arctan — — 

When X t > 0 , added or subtracted 71 when X t < 0 
depending on the sign of Y t , and + 71 / 2. if X t = 0, + Y j . 

Proof. Let ( X t , Y t ) be a solution of the Hopf equations 
( 3 . 7 ). We apply fractional ItO formula to the 
variable R t = -J X 2 + Y~ .Then we have 
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dR > = j[ aX > ~ X ' Y < + £X ‘ i X ' +Y ') + aY ‘ 


+ X,Y t 


+sY, 


T t (4 + 4 )] dt + ^t-[y, 2 {^X, - a 2 Y, f Ht 2H ~ l 


R, 


R: 


+ X; (, * 2 Y, + aft, y Ht 2H ~ { dt -—-\_2Xft (aft, - 


R 

R, k 


-a 2 Y t )(a 2 Y t+ aft t )Ht 2H - l ]dt + ^—[X t (aft t -a 2 Y t ) 
+Y,(a 2 Y t + aft t )]dB? 

- ( aR 2 + sR, 4 ) + -Xr a 2 R, 4+2k - 2 Ht 2H - 1 

R, y ’ R 2 


dt 


1 n 2+fe-l fryH 

+ — a.R, dB, 

R 


= [aR,+ a 2 R, 2k ~ l Ht 2H ~ l + sR 2 ~\dt + cx i R,'dB H 


The argument R, is given by 6 t = arctan 


X,j 


plus 


possibly constants depending on the sign of X t and Y t . 
Applying fractional Ito formula again yields 

dO t = (l - 2(7 ]( 7 2 R; k 2 Ht 2H ' )dt + a 2 R, k \lB H 

which ends the proof of the lemma. 

It is important to note that the equation on the modulus is 
uncoupled of the phase equation on R t The modulus of 

(X„Y.) is therefore solution of a stochastic differential 
equation of type stochastic pitchfork bifurcation . Moreover, 
when k = 1 ,the equations take the simpler form: 

dR - ([a + cr 2 2 Ht 2H ~ l ] R, + sR^dt + aft dB f* 
dO, = (l - 2 a l( j 2 Rt k - 2 Ht 2H - l )dt + cr 2 dB t H 

and hence the variable R t is solution of a stochastic pitchfork 

bifurcation. By application of theorems 3. 4-3.6, it is easy to 
establish the following: 

Theorem 3.7 For k = 1 and: 

(i) H = 1/2 , the null solution of the supercritical Hopf 

equations is almost surely exponentially stable if 

2 2 

a < Ra — Rl and asymptotically stochastically unstable 

2 

cr - -f- cr 2 

if a > — — . In that case, there exists a new 

2 

stochastically stable stationary solution with distribution: 


P a (x) = N a x 


-21 1-4 


V x2/CTi L 


x>0 • 


The null solution of the subcritical Hopf equations is almost 

2 _j_ 2 

surely exponentially unstable if a > — — .It is 


asymptotically stochastically unstable if a > — — + — and 

2 

2 _j_ 2 

stochastically stable if a < l - 2 

2 

(ii) 0 / 1 "C 1/2, the null solution of the supercritical 

Hopf equations is almost surely exponentially stable 
if OC < 0 and asymptotically stochastically unstable if a > 0 . 

(iii) 1/2 < H < 1, the null solution of the supercritical 
Hopf equations is almost surely exponentially stable 
if cr, 2 > cr 2 and asymptotically stochastically unstable 
if a 2 < a\ . 

Remark. Let us emphasize the fact that the case 
0<H <0.5 is substantially different from the general case: 
indeed, we had seen that the singular point was always 
stochastically stable. Surprisingly here the stability of zero 
depends on the relative values of the real and imaginary parts 
of the noise. 

Theorem 3.8 For k < 1 .hi the case of the supercritical 
(s = — l) stochastic Hopf bifurcation, we have: 

(1) H = 1/2 , the null solution of the supercritical Hopf 
equations is almost surely exponentially stable if |cr, | < |cr 2 1 
and asymptotically stochastically unstable if|<X,| > |fTj| . 

(ii) 0 <H <1/2, the null solution of the supercritical 
Hopf equations is almost surely exponentially stable 

if OC < 0 and asymptotically stochastically unstable if a > 0. 

(iii) 1/2 < H < 1, the null solution of the supercritical 
Hopf equations is almost surely exponentially stable 

whatever the parameters cc £ 3? . 

Proof. (1) The case 0 < H < 0.5, falls in the general 
analysis developed in Theorem 3. 5. The general analysis 
hence applies and directly leads to the conclusion of the 
proposition. 

(2) For 1/2 < H <1, Let us consider V(x) — x k as a 
Lyapunov functional for the dynamics. This function is 
clearly C/()R) , i.e. it is twice continuously differentiable 
except X = 0 . Moreover, we have: 

LV< kx k ~' {[a + a 2 Ht 2H ' ] jc - x 3 ) 

-k(l-k)afHt 2H1 x ik - 2 . 

Since we assume that k < 1 and — kx ' ~ at 0 ,the leading 
term close to 0 is of order —k ( 1 — A: ) oy Ht - H 1 x 3k ” which 
is strictly negative for suciently small X whatever the 
parameters OC E $R. 

(3)For H = 1/2. Let us consider V(x) — x k is C 2 ('.K) 
around zero, positive, diverges at zero, we have: 
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LV < kx 


= k 


k - 1 


a + - 


x-x 3 -&(!-&) — x 3i 


( _ 2 _ _ 2 '\ 

<T 2 CTj 


x 3 *- 2 + 


y 


0(x 3k - 2 ). 


and hence is negative closed of zero. 


Theorem 3.9 For k > 1 . In the case of the supercritical 
( s — — l) stochastic Hopf bifurcation, we have: 

(l)When H =1/2, the null solution of the supercritical 
Hopf equations is almost surely exponentially stable 
if a + <J 2 < 0 and asymptotically stochastically unstable 

if a + O' 2 > 0 .In that case, there exists a new stochastically 
stable stationary solution with distribution. 

(2) When 0 < H < 1/2 , the null solution of the 
supercritical Hopf equations is almost surely exponentially 
stable if OC < 0 and asymptotically stochastically unstable 

if a > 0 . 


(3) When 1/2<H<1, the null solution of the 
supercritical Hopf equations is asymptotically stochastically 
unstable whatever the parameters OC e . 

Proof. ( 1 ) When H = 1 / 2 . Let us consider V (x) = x 2 as a 
Lyapunov functional for the dynamics. This function is 
clearly C° ( s -R) ,i.e. it is twice continuously differentiable 


except at X = 0 . Moreover, we have: 


LV <2x 


(a + a 2 ) 


x-x 



Since we assume that k > 1 and — 2x 4 at 0, the term L V 
is equivalent to 2 (a + cr 2 )v 2 and hence locally has the sign 

of a + <7 2 . For a + cr 2 < 0 , we can directly apply 

theorem of [21], For Ct + cr\ > 0 , we use V(x )=— lcg(x) 
again and conclude that 0 is not stable in probability. 

(2) When 0 < H < 0.5 . Let us consider V(x) = X 2 as a 
Lyapunov functional for the dynamics. This function is 
clearly C 2 (]) ?) , i.e. it is twice continuously differentiable 
except at X = 0 . Moreover, we have: 


2x 


\a + a 2 2 Ht 2H ~ l ) 


x-x 


+ 2 a; Ht 


2H-\ 2k 

A . 


^ 21 1 1 21c 

Since we assume that k > 1 and 2 (J"Ht X" at 0, the 


term L V is equivalent to 2 OCX and hence locally has the 
sign of OC . For OC < 0 , we can directly apply theorem 
V.4.1 of [17], X 0 is stable in probability. For (X > 0 , we 

use V(x) = — log(x) again and conclude that X 0 is not 
stable in probability. 

(3) When 0.5 < H < 1 . Let us consider V{x :) = X" as a 
Lyapunov functional for the dynamics. This function is 
clearly C 2 (iR ) , i.e. it is twice continuously differentiable 


except at A" = 0 . Moreover, we have: 


2x 


(a + cr 1 2 Ht 2U 1 ) 


x-x 


. 2 j t . 2H — 1 2k 

+ 2<j x Ht x . 
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Since we assume that k > 1 and — 2x 4 at 0, the leading 
term close to 0 is of order 2 <J 2 HV H l x k which is strictly 
positive for suffciently small X whatever the parameters 

ae'Jl. 


4 Example 

Example 1. Let us consider the Hopfield neural networks 
model[23], classically describing the behavior of a cortical 
column. In that model, the voltage X of a typical neuron in 
the column is solution of the equation: 

dx t = (-ax t + bS (x, ))c// + g(x t ) dB" , (4.1) 

where S(x) is a smooth sigmoidal function describing the 
voltage-to-rate transformation, and b is the typical 
connectivity strength. We choose here S(x ) = tanh( ix) with 
y > 0 the sharpness coeffcient of the sigmoid and 

g(x,) = cr|x| A with k>—. This equation can be locally 

reduced to the normal form of the pitchfork bifurcation since 
the ow is symmetrical, and we obtain using the same 
transformation as performed above: 

drj, = ((- a + b/)r/, + ^ fo . b > y))dt 

+ 2^ k o\rj\ dB" . (4.2) 

When 0.5 1 1 1 , let us denote by Xq the stationaiy 


solution TJ t = 0 a.s. for all t . We have: 

• For k < 1 Tj t = 0 is asymptotically stable in 
probability whatever a,b, y e and (7 X 0. 

• For k > 1 Tj t = 0 is unstable in probability for 
a,b,y e S .K and cr ^ 0 

Example 2. Let us consider a Wilson and Cowan neural 
network[23] composed of an excitatory and an inhibitory 
population. The network equations in that case read: 


dx j = (-Xj (f) + M (x,)-S (x 2 )^df 
+ (cTjXj - cr 2 x 2 ) |x f f * dB , 
dx 2 = (-x 2 (f ) + bS (xj ) + S (x 2 )) dt 
+ (cTjX, - CTjXj ) |x, \ k 1 dB , 


(4.3) 


where S(x) = tanh(/x) with ;/ > 0 and x t =Jxf +x 2 .In 
the deterministic case, it is straightforward to show that the 
system undergoes a Hopf bifurcation at y = 1 . Let us now 

consider the equation satisfied by the modulus p t = | X t | . 
Using fractional It A o's formula we obtain: 


d P,= {- 1+ r)+Y p ’ k l+F ( x 


dt 


+ o lP k dB» 


d0, = 


y + o l a 2 p 2k - 2 + 


O ( p 2 )1 dt + <J 2 p k ~'dB" , 

(4.4) 


where 
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= ‘ + 4) - 2 x 2 x ?) + °(p‘) 

<-f/+0( P ‘). 

When 0.5 < H < 1 , then the dynamics of the system can 
hence be analyzed using the above analysis: 

• For k = 1, the null solution of the supercritical Hopf 
equations is almost surely exponentially stable 


N 

>K 


C |cr 2 . 


• For k < 1, the null solution of the supercritical Hopf 
equations is almost surely exponentially stable whatever 
the parameters y e < '.H. 

• For k < 1, the null solution of the supercritical Hopf 
equations is asymptotically stochastically unstable 
whatever the parameters y e 
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